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ABSTRACT 



> 



Aims. Numerous spectroscopic observations provide compelling evidence for a non-canonical mixing process that modifies the 
surface abundances of Li, C and N of low-mass red giants when they reach the bump in the luminosity function. Eggleton et 
al. (2006) have proposed that a molecular weight inversion created by the ''He(^He,2p)*He reaction may be at the origin of this 
mixing, and relate it to the Rayleigh- Taylor instability. In the present work we argue that one is actually dealing with a double 
^— J ' diffusive instability referred to as thermohaline convection and we discuss its influence on the red giant branch. 

, Methods. We compute stellar models of various initial metallicities using the prescription by Ulrich (1972) (extended to the case 
I of a non-perfect gas) for the turbulent diffusivity produced by the thermohaline instability in stellar radiation zones. 
' Results. Thermohaline mixing simultaneously accounts for the observed behaviour of the carbon isotopic ratio and of the 
t , abundances of Li, C and N on the upper part of the red giant branch. It reduces significantly the '^He production with respect 
' to canonical evolution models as required by measurements of "^He/H in galactic HII regions. 

I Conclusions. Thermohaline mixing is a fundamental physical process that must be included in stellar evolution modeling. 
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d , 1. Introduction ical stellar evolution thcorjoon the RGB. However, spec- 

^ ■ troscopic observations clearly point out that some non- 

• '-H , During the first dredge-up {Idup; Iben 1965), the surface canonical mixing connects the convectivc envelope with 

' composition of low-mass red giant star^ is modified due the HBS and further modifies the surface composition 

. to dilution of the external convectivc stellar layers with of low-mass giants as soon as they reach the bump in 

■ ■ ' hydrogen-processed material: The lithium and the carbon the luminosity function. The sudden drop of the carbon 

abundances as well as the carbon isotopic ratio decline, isotopic ratio ^^C/^^C provides the most pertinent clue 

while the helium 3 and nitrogen abundances increase. The of this mechanism (Gilroy 1989; Gilroy & Brown 1991; 

amplitude of these variations depends on the stellar mass Charbonnel 1994; Charbonnel et al. 1998; Gratton et al. 

and metallicity. This picture is nicely supported by obser- 2000; Shetrone 2003; Recio-Blanco & De Laverny 2007) 

vations. that also modifies the abundances of lithium, carbon and 

nitrogen (e.g., Gratton et al. 2000). This non-canonical 

After the idup the stellar convectivc envelope retreats ' 4- v. ■ ^ -j- a- ^ ^ 1 j. ncor 

,. . . . process appears to be universal as it attects at least 95% 

(in mass) ahead of the advancing hydrogen-burning shell r ^ . u .1 1 -1- c u 4- 

) ' to J to ol the low-mass stars, whether they belong to the held, to 

(HBS) that surrounds the degenerate helium core, and no i 1 1 ^ /m u 1 c t-v m • ^ 

).,„,, 111 open, or globular clusters (Gharbonnel & Do Nascimcnto 

further surface abundance variation is predicted by canon- mnox rri,- u- i, u -t- c -n, i j-- 

^ 1998). ihis high number satisfies the galactic require- 
ments for the evolution of the '^Hc abundance (Tosi 1998; 
Palla et al. 2000; Romano et al. 2003), since the mecha- 
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^ i.e., stars with initial masses below ~ 2 - 2. 5M0 that evolve ^ By this we refer to the modelling of non-rotating, non- 

along the Red Giant Branch (RGB) to high luminosities until magnetic stars, where convection is the only transport process 

helium is ignited in their core under degenerate conditions considered. 
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nism which is responsible for the low values of ^^C/^^C is 
also expected to lead to the destruction of '^He by a large 
factor in the bulk of the envelope material as initialy sug- 
gested by Rood et al. (1984; sec also Charbonnel 1995; 
Hogan 1995; Weiss et al. 1996). 

At the present time, we have no firm physical model for 
the non-canonical RGB mixing. Parametrized approaches 
have been used to reproduce individual observations with 
the diffusive velocity treated as a free parameter (i.e. 
Boothroyd & Sackmann 1999; Weiss et al. 2000). On 
the other hand, rotation-induced mixing has been investi- 
gated thoroughly (see references in Charbonnel & Palacios 
2004). It turns out however that meridional circulation 
and shear turbulence alone do not produce enough mix- 
ing to account for the surface abundance variations as re- 
quired by the observations (Palacios et al. 2006). 



2. Mixing due to thermohaline convection 

Recently Eggleton et al. (2006, 2007) identified a possible 
cause for such non-canonical mixing, namely the molecular 
weight inversion created by the •^He('^He,2p)"'He reaction 
in the external wing of the HBS. As it was pointed out 
already by Ulrich (1972), this nuclear reaction is somewhat 
singular in that it produces more particles per unit mass 
than it started from. Using their 3D hydrodynamic code 
to model a low-mass star at the RGB tip (Dearborn et 
al. 2006), Eggleton and collaborators found that such a 
/i-profile leads to efficient mixing, of the kind required to 
"reconcile Big Bang and stellar nucleosynthesis (as far as 
^He is concerned)", as they put it. 

According to them, the instability responsible for that 
mixing is the well-known Raylcigh- Taylor instability, a dy- 
namical instability which is triggered when a layer of heav- 
ier fluid lies over lighter fluid. In stellar interiors, that in- 
stability takes the form of the convective instability, which 
tends to render the temperature gradient adiabatic rather 
than to suppress the density inversion. But what first oc- 
curs in a star, as the inverse /i-gradient builds up, is actu- 
ally a double diffusive instability, which was discussed long 
ago in the literature under the generic name of "thermo- 
haline convection" (Stern 1960). This instability appears 
in various astrophysical situations, for instance when ^He 
or C-rich material is deposited at the surface of a star 
in a mass transferring binary (Stothers & Simon 1969; 
Stancliffe et al. 2007). Recently it has been invoked when 
a star accretes heavy elements during planet formation 
(Vauclair 2004). 

The thermohaline instability differs from the convec- 
tive instability in that it involves two components, of 
which one, the stabilizing one (temperature) diffuses faster 
than the other (salt), whose stratification is unstable. It 
occurs in a stable stratification that satisfies the Ledoux 
criterion for convective stability: 



but where the molecular weight decreases with depth: 



'ad 



V 



dln/i 
dlnP 



< 0. 



(2) 



We use the classical notations for V = {d\nT/dlnP), 
ip = {d\np/d\nfi)p^T and S = — (91np/c)lnr)p^^. 

In the laboratory, the instability takes the form of "salt 
fingers"; since heat diffuses faster than salt, these fingers 
sink because they grow increasingly heavier than their en- 
vironment, until they become turbulent and dissolve. In 
stellar interiors, the role of salt is played by a heavier 
species, such as helium, in a hydrogen-rich medium. 

Ulrich (1972) was the first to derive a prescription for 
the turbulent diffusivity produced by that instability in 
stellar radiation zones. Through a linear analysis, and as- 
suming perfect gas {ip ^ S = 1) he gets (his eqs. 24 and 
31) 



Dt = Ct K 



(3) 



K is the thermal diffusivity. His non-dimensional coeffi- 
cient involves the aspect ratio a (length/ width) of the 
fingers: 



(4) 



for the value he advocates, a — 5, this coefficient is rather 
large: Ct = 658. 

Kippenhahn et al. (1980) extended Ulrich's expression 
to the case of a non-perfect gas (including radiation pres- 
sure, degeneracy): 



Dt = CtK 



(I) 



6 J (Vad - V) 



for V„ < 0. 
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V,, > 0, 



(1) 



However instead of fingers they considered "blobs" which 
according to them are destroyed very rapidly: their mixing 
length is of the order of the blob size and this translates 
into a smaller value of 12 for the coefficient Ct- 

It is of course disturbing that these coefficients Ct differ 
by almost two orders of magnitude, depending on whom 
one believes; it would be most desirable to confront them 
with realistic 3D simulations. In the meanwhile, we tend 
to favour a large value for the coefficient Ct as advocated 
by Ulrich, because all experiments so far have displayed 
slender fingers, rather than blobs (cf. Krishnamurti 2003). 
Crude as it may be, this prescription has the advantage of 
being rooted in the physical process, unlike the admittedly 
ad hoc diffusivity chosen by Eggleton et al. (2007). 

3. IVIodel calculations 

We compute three evolution models of a 0.9 Mq star with 
various initial values of [Fc/H], namely -1.8, -1.3 and -0.5, 
with the aim of comparing our theoretical predictions with 
observations over a large range of metallicity. We use the 
code STAREVOL (Siess et al. 2000; Palacios et al. 2003, 
2006). In the present study the transport of particles in 
the radiative regions is due to thermohalinc mixing only. 
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Fig. 1. Profiles of the logarithm of the mass fraction of H, 
^He, i^C, i^C, "N, 1^0, and of the mean molecular weight 
gradient = d In^/ d InP (bold line, full when > 
0, dashed otherwise) as a function of the reduced mass 
(see text) inside a 0.9 Mq, [Fe/H]=-1.3 star. The left and 
right panels correspond to models including thermohaline 
mixing and located just before and after the bump. In 
the latter case, the thermohaline instability has just set in 
naturally, while in the former case it is still prevented due 
to the strong /i-gradient left behind by the Idup 

The expression we use for Dt is that given in Eq. [5] that 
includes the correction for a non-perfect gas. For simplicity 
we assume Ct = 1000, which is of the order of magnitude 
of the coefficient advocated by Ulrich (1972). The model at 
[Fe/H]=-1.3 is also computed with Ct = 100 and 12. We 
do not consider any additional transport process related 
to rotation. 

Fig. [T] presents the abundance profiles of selected el- 
ements and the mean molecular weight gradient = 
d In /J,/ dlnP, just before and after the bump in the 0.9 
M0, [Fe/H]=-1.3 star when Ct=1000. The abscissa is the 
relative mass SM ~ {Mr — A/coic)/(^^cnv ~ -Aicorc) defined 
as ranging from to 1 between the bottom of the HBS 
and the base of the convective envelope. 

During the idup phase, the convective envelope ho- 
mogenizes the star down to very deep regions, and builds 
a very steep gradient of molecular weight at the point of 
its maximum penetration. This corresponds to the exter- 
nal peak (at SM ~ 0.12) in the profile of that can 
be seen in the left panel of Fig. [T] On the other hand 
the deeper peak (at SM ^ 0.05) locates the region where 
H is efficiently depleted by nuclear reactions in the HBS. 
One sees that at the current luminositity (L=9OL0) both 
peaks are very close to each other, meaning that the star 
is approaching the bump. At this evolutionary point, 



is positive in the whole radiative region, and the thermo- 
haline instability cannot set in (see Eq. 2). The profiles of 
the chemical elements are thus identical to those obtained 
in a canonical model. 

When the HBS passes through the /Lt-discontinuity left 
behind by the Idup (i.e., at the bump), H-burning occurs 
in a homogeneous region. The ■^He('^He,2p)"'He reaction 
slightly lowers the molecular weight in the external wing 
of the HBS where becomes negative (bold dashed line 
in the right panel of Fig. [T|), allowing the thermohaline 
instability to develop naturally between the ^He-burning 
region and the base of the convective envelope. Deeper 
inside the radiative region, remains positive (bold full 
line) and no thermohaline mixing occurs. As can be seen 
in Fig. [2 the surface abundance of ^He, ^^C, i^C, and 
are already modified soon after the onset of thermohaline 
mixing. However the surface abundance of ^^O remains 
constant because the thermohaline mixing does not extend 
down to the very deep region where full CNO-burning 
operates at equilibrium. 

We show in Fig.[2]the profile of the diffusion coefficient 
Dt (Eq. 5, for Ct=1000), for the same model just after 
the onset of the thermohaline instability (L=115L0). Dt 
is high enough to connect the HBS and the convective 
envelope and to modify the surface abundances as shown 
below. It is two to three orders of magnitude higher than 
that characterizing the rotational mixing (Palacios et al. 
2006). 

4. Signatures of mixing and comparison with 
observations 

In Fig. [3] we compare the temporal evolution of ^Li, 
i^C/i^C, [C/Fe], and [N/Fe], obtained at the surface of 
our models with homogeneous observational data for field 
stars with -2<[Fe/H]<-l by Gratton et al. (2000). We see 
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Log (L/Lg) 

Fig. 3. Evolution of the lithium abundance, of the car- 
bon isotopic ratio ^^c/i^C, of [C/Fe] and [N/Fe] as a 
function of the luminosity logarithm for the O.OM© mod- 
els with [Fe/H]=-1.8, -1.3 and -0.5 (respectively dashed- 
dotted blue, full line black, and dashed red) computed 
with Ct = 1000 in Eq. 5. The black dotted line is for the 
model with [Fe/H]==-1.3 calculated with Ct=100. The ar- 
rows in the upper panel indicate the location of the bump 
for the three metallicities. Observational data are from 
Gratton et al. (2000) for field stars in the metallicity 
range [Fe/H] G [-2;-l]. Blue symbols are for stars with 
measured [Fe/H] values lower than -1.4, black for stars 
with -1.4 <[Fe/H]< -1.2, and red for [Fe/H] higher than 
-1.2. Circles are actual measurements, open upward tri- 
angles are lower limits and open downward triangles are 
upper limits. 

first that variations associated with the Idup fit well the 
observational behaviour on the lower RGB. The surface 
abundances then stay constant until the stars reach the 
bump. Until that evolutionary point the predictions arc 
identical to those of the canonical models. 

Then at the luminosity of the bump (which increases 
for a given stellar mass when the initial metallicity de- 
creases, as shown by the arrows) thermohaline mixing 
leads to a second episode of abundance variations that ex- 
plains remarkably well the datgffl Li is rapidly destroyed, 
and the carbon isotopic ratio reaches values between 5 
and 7 very close from the equilibrium value. [C/Fe] and 

^ Oscillations in stellar luminosity can be seen just above the 
bump. They are related to the transport of ^He from the con- 
vective envelope into the outer HBS region where this element 
burns in favour of ^Be. The related nuclear energy release is 
large. When it happens to rival that of the pp-chains, the total 
luminosity increases and the star readjusts. Once the abun- 
dance of ^He is low enough the oscillations vanish. 



[N/Fe] respectively decrease and increase. These C and N 
variations are stronger for lower initial metallicity, as re- 
quired by the observations. As already noted in § 3, the 
surface abundances of the heavier elements (^^O and ^'^Na 
in particular, not shown here) which arc affected by nu- 
clear reactions much closer to the He-core do not vary at 
the surface. 

The models shown in Fig. [3] were computed assuming 
Ct = 1000 in Eq. 5. We also show the predictions for the 
model with [Fc/H]=-1.3 computed with Ct = 100. In that 
case the surface abundances change less rapidly and one 
does not reproduce the abrupt variation at the luminosity 
of the bump depicted by the observations. Our analysis is 
thus compatible with the high value for the coefficient Ct 
advocated by Ulrich (1972) and it excludes the low value 
(12, not shown here) derived by Kippenhahn ct al. (1980). 

5. The case of '^He 

On the main sequence, a "^He peak builds up due to pp- 
reactions inside the low-mass stars (Ibcn 1967), and is 
engulfed in the stellar envelope during the Idup. As a con- 
sequence the surface abundance of "^Hc strongly increases 
on the lower RGB as can be seen in Fig. [4] for our models 
at three different metalicitics. Its value reaches a maxi- 
mum when the whole peak is engulfed. After the Idup 
the temperature at the base of the convcctivc envelope is 
too low for ■^He to be nuclearly processed. As a result in 
canonical models ■^He stays constant at the surface and its 
final value is strongly increased with respect to the initial 
one (this is the value before thermohaline mixing sets in 
at the bump). 

After the bump, thermohaline mixing brings '^He from 
the convective envelope down to the region where it is nu- 
clearly burned. This leads to a rapid decreases of the sur- 
face abundance of this clement as can be seen in Fig. 21 
This confirms the early suggestion by Rood et al. (1984) 
that the variations of the carbon isotopic ratio and of 
■^Hc are strongly connected (see also Charbonnel 1995; 
Egglcton ct al. 2007). It is important to note that in the 
models presented here "^He decreases by a large factor in 
the ejected material with respect to the canonical evolu- 
tion predictions but that low-mass stars remain net pro- 
ducers (while far much less efficient than in the canon- 
ical case) of '^He. Computations for different masses and 
metallicities have now to be performed in order to estimate 
the actual contribution of low-mass stars to Galactic ■^Hc 
in the framework proposed here. We are confident that 
the corresponding '^He yields will help reconciling the pri- 
mordial nucleosynthesis with measurements of ''He/H in 
galactic HII regions (Balser et al. 1994; Bania et al. 2002; 
Charbonnel 2002). 

6. Conclusions 

It is known since a long time that a non-canonical mixing 
process modifies the surface abundances of low- mass stars 
when they arrive at the bump in the luminosity funtion. 
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Fig. 4. Evolution of the surface abundance of ^He (in mass 
fraction). The hnc symbols are as in Fig.[3l 

Based on SD-modcling of a low-mass star at the RGB tip, 
Eggleton and collaborators have proposed that the molec- 
ular inversion created by the '^He('^He,2p)*He reaction in 
the external wing of the HBS may be the root cause for 
such a mixing. They ascribe this mixing to the Rayleigh- 
Taylor instability. 

In the present work we argue that one is actually deal- 
ing with a double diffusive instability called thermohaline 
convection, which has been discussed long ago in the liter- 
ature. We compute stellar models including the prescrip- 
tion given by Ulrich (1972) and extended to the case of 
a non-perfect gas for the turbulent diffusivity produced 
by that instability in stellar radiative zone. The results 
presented here indicate that thermohaline convection si- 
multaneously accounts for the observed behaviour of the 
carbon isotopic ratio and for the abundances of Li, C and 
N in RGB stars at and above the bump in the luminos- 
ity function. It also avoids large "^He production by low- 
mass stars as required by chemical evolution models of 
the Galaxy. It does not modify the O nor the Na surface 
abundances. In a future work the effect of the thermoha- 
line mixing in models of various masses will be presented. 
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